A multiple-scattering approach to interatomic interactions and superradiance in 

inhomogeneous dielectrics 
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The dynamics of a coUection of resonant atoms embedded inside an inhomogeneous nondispersive 
and lossless dielectric is described with a dipole Hamiltonian that is based on a canonical quantiza- 
tion theory. The dielectric is described macroscopically by a position-dependent dielectric function 
and the atoms as microscopic harmonic oscillators. We identify and discuss the role of several 
types of Green tensors that describe the spatio-temporal propagation of field operators. After in- 
tegrating out the atomic degrees of freedom, a multiple-scattering formalism emerges in which an 
exact Lippmann-Schwinger equation for the electric field operator plays a central role. The equa- 
tion describes atoms as point sources and point scatterers for light. First, single-atom properties 
are calculated such as position-dependent spontaneous-emission rates as well as differential cross 
sections for elastic scattering and for resonance fluorescence. Secondly, multi-atom processes are 
studied. It is shown that the medium modifies both the resonant and the static parts of the dipole- 
dipole interactions. These interatomic interactions may cause the atoms to scatter and emit light 
cooperatively. Unlike in free space, differences in position-dependent emission rates and radiative 
line shifts infiuence cooperative decay in the dielectric. As a generic example, it is shown that near 
a partially reflecting plane there is a sharp transition from two-atom superradiance to single-atom 
emission as the atomic positions are varied. 

PACS numbers: 42.50.-p, 41.20.Jb, 42.70.Qs 



I. INTRODUCTION 

The spontaneous-emission rate of an atom depends on 
its dielectric environment P, Q] and in particular on the 
precise position of the atom if the medium is inhomoge- 
neous IE 0) IE IE 01 ■ Spontaneous emission can only 
be understood quantum mechanically, but the classical 
Green function determines the emission rate. In partic- 
ular, the emission rate is proportional to the imaginary 
part of the Green tensor of the medium at the atomic 
position. The dipole-angle average of the emission rate 
is also known as the local optical density of states 0] . 

In order to study the influence of the medium on much 
more than just single-atom spontaneous-emission rates, 
in this paper a rather general multiple-scattering theory 
is set up. It is is based on 'macroscopic quantization' 
theories of the electromagnetic field in inhomogeneous 
lossless dielectrics, see for example fs", (IE IH El ■ Mi- 
croscopic treatments of (macroscopically homogeneous) 
dielectrics in quantum electrodynamics can be found in 
|l3lll4j but will not be used here. An important result in 
the macroscopic theories is that photons can be defined 
as the elementary excitations of the true modes of the 
dielectric. Corresponding mode functions are the (classi- 
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cal) harmonic solutions of the wave equation. 

Emission rates of an atom not only change due to the 
nonresonant dielectric environment, but also due to the 
presence of other atoms with the same transition fre- 
quency. As is known since the pioneering work by Dicke 
[la], resonant atoms in each others neighborhood decay 
cooperatively. Depending on the many-atom state, the 
atoms decay faster than a single atom up till twice the 
single-atom rate (superradiance) or decay slower or not 
at all (subradiance) . Lifetime changes of individual atom 
pairs as a function of their distance were measured only 
recently IT^ ; for two Ba^ ions that emit at a wavelength 
of 493 nm and for well-defined separations |R| around 
1.5 /im, subradiant and superradiant lifetime effects of 
less than ±2% were observed. 

Superradiance occurs for the so-called Dicke states that 
have a zero expectation value of the total dipole moment 
EEli but also for atomic product states with a nonzero 
dipole moment 0, 0|. Superradiance also occurs for 
classical dipoles. It is a general phenomenon also exhib- 
ited in acoustics by nearby identical tuning forks, or by 
strings in a piano ^El- These systems have in common 
that the atoms (or oscillators) interact with a field that is 
influenced by the radiation reactions of all nearby atoms 
together. 

Cooperative effects of resonant atoms will be influ- 
enced by their dielectric environment. In this paper the 
influence of a nondispersive and lossless inhomogeneous 
dielectric on embedded or nearby resonant atoms is stud- 



2 



ied. Quantum theory is used both to describe the Hght 
and the atoms. As we focus on the effects of the inho- 
mogeneous dielectric, the atoms are modelled simply as 
quantum harmonic oscillators in their ground states or 
first excited states, with fixed dipole orientations. To 
be sure, in choosing this model we neglect optical satu- 
ration effects of the atoms. The dielectric is described 
macroscopically in terms of a real-valued relative dielec- 
tric function £(r), the form of which will be left arbitrary. 
The precise measurement of two-atom superradiance in 
free space as a function of distance is a fundamental test 
for quantum electrodynamics ■ The effects calculated 
here are a test for macroscopic quantization theories for 
inhomogeneous dielectrics. 

Like single-atom emission [T9l |. two-atom superradi- 
ance will be modified in the close vicinity of a mirror, 
or inside an optical cavity 20]. Strong modifications of 
superradiance are also predicted for photonic crystals, di- 
electrics with periodic refractive-index variations on the 
scale of the wavelength of light 0, |2^ . Large effects 
are predicted in the so-called isotropic model for a pho- 
tonic crystal, which is really a toy model in the sense that 
all local an orientational inhomogeneities of the electro- 
magnetic field are neglected. In a real photonic crystal, 
two-atom superradiance is expected to sensitively depend 
on the coordinates of both atoms. The present formal- 
ism is valid for an arbitrary real dielectric function and 
encompasses the interesting special cases just mentioned. 

Atoms that exhibit superradiance interact strongly 
enough to share and exchange the optical excitation be- 
fore emission. The more common and better studied 
situation for resonant atoms in a dielectric is that the 
interaction between the atoms is weak compared to in- 
teractions with baths that the individual atoms have. 
Then the optical excitations are transferred irreversibly 
from donor to acceptor atoms via a process called 'res- 
onance energy transfer', as described by Forster theory 
[2^ and its modern generalizations Resonance en- 
ergy transfer is influenced by the dielectric environment. 
For example, calculations show that two-atom interac- 
tions can be strongly influenced by an optical microcavity 
HE nil I since the cavity modes with eigenfrequencies 
close or equal to the atomic transition frequency play 
a dominant role. Indeed, experiments have shown that 
the interatomic (dipole-dipole) interaction is increased 
when the atoms are placed in a cavity at positions where 
resonant optical modes have their maxima (2^]. In an- 
other interesting experiment, Forster excitation transfer 
is found to scale linearly with the local optical density 
of states at the donor position 30] . Although we focus 
on superradiance, the formalism in the present work is 
quite general and can also be used as a quantum elec- 
trodynamical foundation for the study of energy transfer 
processes in inhomogeneous media. Recent progress in 
this direction can be found in |^ [s^ . 

Often in quantum optics an 'all-matter' picture is em- 
ployed, where the dynamics of the electromagnetic field 
is integrated out, for example in the optical Bloch equa- 



tions |17l l33| . Here instead we treat spontaneous emis- 
sion and superradiance in an 'all light' picture, which is 
convenient when studying the effect of the dielectric. A 
multiple-scattering theory is set up in which the atoms 
show up both as sources and as scatterers of light. It 
is known that superradiance can be viewed as caused by 
multiple-scattering interactions |3 Issl |3p] . Light scat- 
tered off a collection of atoms will show multi-atom reso- 
nances and cooperative effects, also due to multiple scat- 
tering. In quantum scattering theory such resonances 
appear as well and sometimes are called 'proximity reso- 
nances' mni]. 

The concept of a point scatterer proved very fruitful in 
the study of multiple-scattering of classical light in free 
space [33, |4l| . Here, the point-scattering formahsm 
will be put to use in quantum optics of inhomogeneous 
dielectrics. Multiple light scattering will be described in 
terms of Green functions of the medium. The emphasis 
of the paper will be on formalism, but it ends with an 
application to superradiance in a model dielectric. 

The paper has the following structure: in Sec. [H] the 
point-scattering model for interacting guest atoms is in- 
troduced. Properties of several types of Green functions 
of the medium are discussed in Sec. lIIII Sec. IIVI discusses 
medium-induced modifications of single-atom properties 
such as spontaneous-emission rates and elastic scattering. 
The generalization to a finite number of host atoms is dis- 
cussed in Sec. The formalism is applied to two-atom 
superradiance in Sec. IVII in particular to superradiance 
near a partially reflecting plane in Sec. IVIII Conclusions 
are drawn in Sec. IVIIII 



II. ATOMS AS POINT SOURCES AND AS 
POINT SCATTERERS 

A. The Hamiltonian 

Consider an inhomogeneous dielectric with relative di- 
electric function e(r) with a finite number of embedded 
neutral atoms. The dipole Hamiltonian for this system 
is the sum of a field part, an atomic part, and an inter- 
action part between field and atoms. More precisely, the 
Hamiltonian can be found after canonical quantization 
[Hill to have the form H = Hp + Ha + Haf, with 

Hf = ''^fvjj\a\a\ (la) 

A 

N 

Ha = ^ftr!™6|„6™ (lb) 

m— 1 

N 

Uaf = -^/x™-r(R„) 

m— 1 

= ^{bni + h\^){gx„^ax+ gl„,a{). (Ic) 
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Notice that there is no direct interaction term between 
neutral atoms. In a minimal-coupling Hamiltonian, there 
would have been such a direct coupling term. The sit- 
uation is analogous to the free-space case [s^, E3 • The 
field part Tip of the Hamiltonian is a sum (or integral) 
over harmonic oscillators corresponding to the harmonic 
solutions ('true modes') of the Maxwell equations for 
the inhomogeneous dielectric in the absence of the atoms: 

-V X V X fA(r) +e(r)(LJA/c)2fA(r) = 0. (2) 

For oja 7^ 0, these modes are generalized transverse, 
which means that V- [£(r)fA(r)] = 0. Their orthonormal- 
ity condition reads / dre(r)fj^(r) • fA'(r) = 6\\', where * 
denotes complex conjugation. The modes are complete, 
in other words they form a basis for the subspace of gen- 
eralized transverse functions. For free space [e(r) = 1] 
the fA are the well-known transverse plane-wave modes. 

In the atomic Hamiltonian TCa, the atomic transition 
frequencies flm and transition dipole moments fim may 
be all different, either because the guest atoms are of 
different species or because identical atoms feel a differ- 
ent environment. The frequencies Qm are assumed real, 
which means that nonradiative broadening is neglected. 
The atoms are very simply described as harmonic oscilla- 
tors with frequencies This is a good approximation 
within a certain frequency range and as long as satura- 
tion effects of the upper atomic state can be neglected. 
The atomic transition dipole moments fim are assumed 
to be real- valued and to have fixed orientations. This 
assumption is better for molecules or quantum dots in a 
solid surrounding than for atoms in the gas phase. For 
convenience, the name 'atoms' will be used for the guests 
in the dielectric. The operators bl^{t) create atomic exci- 
tations by annihilating an atom in the ground state while 
at the same time creating the atom in the excited state. 

The total displacement field D(r, i) is equal to the 
displacement field eo£(r)E(r, t) of the inhomogeneous 
medium plus the sum J2m'^rn.i''^,t) of the polarization 
fields produced by the guest atoms. In the dipole ap- 
proximation, these polarization fields have the form 

P™(r,t) - d{r-R,n)Pm{t) 

= ^(r - R„) fim [ h,n{t) + blit) ] . (3) 

In the dipole interaction term T^af of the Hamiltonian, 
a field called F was introduced that is an abbreviation of 

F(r,t)^D(r,t)/[eoe(r)]. (4) 

Atomic dipoles couple to this field F(r, i) It is 

equal to the electric field operator E(r, t) everywhere, ex- 
cept at the positions R™ of the guests, since the guest 
dipoles couple to fields in which their own polarization 
fields are included. For free space this self-interaction in 
the dipole coupling is known [4^ . The mode expansion 
of the field F(r, t) has a simple form, being the sum of a 



positive-frequency part F^"*"^ (r, t) containing only annihi- 
lation operators and its Hermitian conjugate F^^^(r,t), 
where 

FW(r,t)=z^ygaA(t)fA(r). (5) 

In the absence of the atoms, the time dependence of 
the annihilation operators in (|SJ) would be harmonic and 
F(r, t) would be equal to the electric field E'^°) (r, t). Here 
and below, the superscript (0) denotes the absence of 
guest atoms in the inhomogeneous dielectric. For conve- 
nience, coupling constants between atom m and optical 
mode A in Eq. l(Tc|) are defined as 



9X771 = -iy ^^fJ'ni ■ hCRm)- (6) 

Notice that the coupling constants g\m are zero for (lon- 
gitudinal) modes corresponding to wa = 0. It is by a 
convenient choice of gauge that the longitudinal modes 
are decoupled from the atoms in the Hamiltonian 



B. Derivation of Lippmann-Schwinger equation 

The goal of this section is to derive a Lippmann- 
Schwinger equation for the field F inside the inhomoge- 
neous dielectric in the presence of the N guest atoms, by 
integrating out the atomic dynamics. Heisenberg's equa- 
tion of motion leads to the following equations of motion 
for the field operators: 

dA = -iivxax- (i/tij^gl^ib^ + bl^) (7a) 

d[ = iujxa{ + {i/n)^gxr7i{b,7i + blJ. (7b) 

rn 

(The dot denotes the time derivative; explicit time de- 
pendence of the operators is henceforth dropped.) The 
field operators are coupled to the atomic operators and 
the operators of atom m satisfy the equations 

b„, = -i^77ibm- {ilK)^{gxma\ + glmAm) (8a) 

A 

bL = i^mb\^ + [i/h) ^(ffAmflA + ^AmflL)- (8b) 
A 

Now take the Laplace transform (or one-sided Fourier 
transform) of the equations of motion. The transform 
will have the argument — ioj, for example bm{^) = 
dt e^^^hmit). Here and in the following the frequency 
Lj is assumed to contain an infinitesimally small positive 
imaginary part so that the transform is well-defined. The 
equations are algebraic after the transformation. 

Also in Fourier language, the equations for the 
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frequency-dependent atomic operators become 

ibrnit = 0) 



^\^gxmax{Lu) + gl^al^Juj) (9a) 



iblit = 0) 



LU + fir 



^\gx,nax{uj) + g*x^a\^{uj) (9b) 



In obtaining these equations, it was assumed that at time 
zero, the annihilation operators ax{t) coincide with the 
a^^\t), the operators in the absence of the guest atoms. 



The latter operators have the simple harmonic time de- 

0, the transform of which 



pendence a^°^ {t) - 
becomes —i{Lj 



-iu!xaf\t) 



(0) 



(c) 



,(0) 



[t = 0) after a partial 



integration. Notice that 6m (w) and h'l^{ijj) in Eq. ^ are 
defined as the transforms of bm{t) and bl^{t), respectively. 
The time-dependent operators are hermitian conjugates 
(6j'„(i) = [bm{t)]^), but the frequency-dependent opera- 
tors are not (fe|„(t^) ^ [b,n{uj)]^). 

The right-hand sides of Eqs. H9al) and (|9b|) will now be 
used to replace bm{(^) and blj^{uj) in the Laplace trans- 
forms of the Eqs. H7a|) and Ij7bp for the field operators. In 
doing this, the atomic dynamics is integrated out. One 
obtains for the frequency-dependent annihilation and cre- 
ation operators of the electromagnetic field 



axiio) 
a{{Lj) 



- JO) 



UJ - LUX 
UJ + UJX 



'y^,9xm 

ni 



bnM , 5t„(0) 



bm{0) , 6t„(0) 



E 



UJ — UJX — , UJ 
m^X 



A m"™ 



E 



7n.,\' 



gx'ra aA'(w) +gvm al'(^) ilOa) 
5A'™ ax'[uj)+g*yra oU'^ljob) 



r 



The optical modes are no longer independent because of 
the interaction with the atoms. The three terms in the 
right-hand sides of Eqs. (|10a|l and Ijl0b|) can be related to 
three reasons why there can be light in mode A: firstly, 
because there is light in the undisturbed mode that has 
not 'seen' the atom; secondly, because the atom can emit 
light into the mode A; the third term describes transitions 
of light in and out of the mode A to and from modes A', 
due to scattering off one of the guest atoms. Since the 
relations (|10|l are implicit rather than explicit solutions 
for the operators, the identification of terms in the equa- 
tions with scattering and emission processes can only be 
approximate. 

The results (|10ll for the creation and annihilation op- 
erators can be directly used with Eq. to find the fol- 
lowing equation for the field F 



F(r,w) = E(°)(r,w) 

+ ^ K(r,Rm, w) • Sm(t^) 



(11a) 
(lib) 



^ K(r, Rm, w) • Vm(c^) • F(Rm, w)(llc) 



This is the central result of this paper. It is an ex- 
act Lippmann-Schwinger equation and it describes the 
resonant scattering off and emission by guest atoms in- 
side an inhomogeneous dielectric, both for strong and for 
weak atom-field interactions. The equation has an undis- 
turbed term (|lla|l . a source term (|llb|l . and a scattering 
term (|llcll . 



The elements of Eq. (|11|) must still be explained. The 
operator E'*'^ (r, uj) is the electric field in the absence of 
the atoms, with both the positive and negative frequency 
parts. The atomic source operators Sm(w) are vectors 
that have the form /im'5'm(w), where fim denotes the unit 
vector in the direction of the atomic dipole moment /Xm 
and 



eoc^ 



bmiO) , bliO) 



(12) 



Notice that Sm features the atomic creation and anni- 
hilation operators at the initial time zero: in quantum 
optics, the atomic variables can not be completely inte- 
grated out in an 'all-light' picture. 

The optical potentials Vm(w) produced by the atoms 
are dyadics equal to fi„iVm{uj) ft-rm where 



Vm{uj) 



7 ? 
Mm^ 

heoc^ 



1^2 



(13) 



Both the sources and the potentials have resonances 
at frequencies ifim- Potentials Vm{uj) are sometimes 
rewritten as — (cj/c)^ times a 'bare polarizability' aBmiuj) 
. In the present case, the bare polarizabilities are real 
(except exactly on resonance) and they change sign when 
going through their resonances at f2m; the resonances are 
infinitely sharp because all possible nonradiative decay 
processes are neglected; the polarizability is called 'bare' 
because it does not (and should not) contain radiative 
broadening of its resonance (but see Sec. II V|). 
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The last undefined factor in Eq. is the dyadic 

quantity K which is given by 



K(r,r',c^)=c' 



A 



fA(r)f,*(rO 



(14) 



Usually, in a Lippmann-Schwinger equation one finds the 
Green function (called G) of a medium where we now find 
the dyadic K. Interestingly, K turns out to be different 
from G, even for free space, as will be studied in Sec. lIIII 
All the elements of Eq. have now been defined. 

Another important field operator for the medium is the 
vector potential A. The magnetic field B equals V x A. 
In the canonical quantization theories [lllll2j | upon which 
our Hamiltonian (QJ is based, the generalized Coulomb 
gauge is chosen, which means that A is generalized trans- 
verse. Its expansion in terms of the normal modes is given 
below. With Eq. (jTUI) this leads to 



A(r,c.) ^ J2 



aA(^)fA(r)+4(^)fi:(iOLL) 



+ T^y'G'^(r,R,„,cj) • S„i(cj) 

m 
m 

Analogously to Eq. an undisturbed term, a source 

term and a scattering term can be identified for the vector 
potential. 

A difference between Eq. (Illll for the field F and 
Eq. l)15bp for A is that only the former is a Lippmann- 
Schwinger equation and that A immediately follows 
from the solution of F, rather than vice versa. In a 
minimal-coupling formalism, one would find a Lippmann- 
Schwinger equation for the vector potential instead. An- 
other important difference between the equations for the 
two fields is that in Eq. (|15b|l for A the generalized trans- 
verse Green function G appears, rather than the dyadic 
K of Eq. Hll|) . Definitions of and relations between G, 
G"^, and K will be studied shortly, in Sec. IIIII 

Often, Lippmann-Schwinger equations are derived in 
'all-light' formalisms that start with a given optical po- 
tential as a perturbation. Here instead, the approach 
started one level deeper and the optical potential 
is output rather than input. An important feature in 
Eq. Hll() is that the atoms are not only point scatterers 
(potentials), but also point sources for light. Both appear 
as two sides of the same coin in one equation. Solutions 
for the equation will be discussed shortly in Sec. IIVI for 
one atom and in Sec. for several atoms. 



III. GREEN FUNCTIONS OF THE MEDIUM 

The dyadic quantities K and G"^ will now be related 
to the Green function of the medium. The (full) Green 



tensor G(r,r',u;) of an inhomogeneous medium charac- 
terized by the dielectric function e(r) is the solution of 
the wave equation 

-V X V X G(r, r', lj) + e{r){uj/c)^G{r, r', oj) = 6{r - r')l, 

(16) 

where the right-hand side is the ordinary Dirac delta 
function times the unit tensor. 

For a discussion of G""", it is useful to first introduce 
the concept of a generalized transverse delta function 
IT^ . (For comparison. Green and delta functions of a 
homogeneous medium are given in the Appendix.) A 
generalized transverse delta function dj (a distribution) 
can be defined in terms of the mode functions f\ [see 
Eq. ©]: 



^fj:(r)f,(r')e(r'). 



(17) 



Now dJ has the projection property / dri dj{ri,r) ■ 
X''"(ri) = X'^(r) for all (ordinary) transverse vector 
fields X""". The bar in 6j denotes the transpose. The 
same projection can be applied to Eq. (|15|) . In doing so, 
the transverse double-curl term is projected onto itself. 
The generalized transverse Green function G"*" can now 
be defined such that e(r)G"''(r, r', a;) equals the projec- 
tion J dri^J(ri,r) • [e(ri)G(ri, r')]. The projection then 
leads to the following equation for G""": 

-Vx VxG(r.r'.tj) + e(r)(w/c)2G'^(r,r',w) = ^J(r',r). 

(18) 

Notice that G rather than G""" appears in the first term. 
Furthermore, a longitudinal Green function G'^ can be 
defined as G — G"^. By taking the difference of Eq. Hlt)|) 
and Eq. (|18l) one can see that G has the form 



G^(r,r'^ 



1 



e{r){uj/c)^ 
1 



[J(r-r')l -^J(r',r)](19a) 



(19b) 



In equality Ijl9b|l the generalized longitudinal delta func- 
tion was defined as the difference between the ordi- 
nary Dirac and the generalized transverse delta function, 
so that dJ + dl' = 6\. We called G'" the longitudinal 
Green function, but it is not self-evident that for ev- 
ery inhomogeneous dielectric G'" is longitudinal indeed. 
Proofs that / dr'G^(r,r') • X'^(r') = and also that 
J drX"'^(r) • G'"(r,r') — can be found with the help 
of Eqs. (32a) and (32b) of Ref. Ull, respectively. Then, 
since G'" is longitudinal, G in Eq. H18|) can be replaced by 
G"^. Hence the projection of Eq. I|16|l leads to a unique 
defining equation for G"^. 

From Eqs. Q and H18(l . it follows that the generalized 
transverse Green tensor G""" has the mode expansion 



fA(r)fj:(rO 

(liJ + ir\)^ — uj'^ 



(20) 



6 



In this manifestly generalized transverse form, G"*" ap- 
peared in Eq. Ijl5b|) for the vector potential. In the de- 
nominator of Eq. (|2()|l we have for once made explicit 
the positive and infinitesimally small imaginary part of 
the frequency oj, through the term irj. With the positive 
sign of the imaginary part, H2()|l is the causal Green func- 
tion which transformed back to the time-domain gives a 
Green function G'^(r, r', t — io) which is nonzero only for 
positive time differences — tp). 

We are now in the position to rewrite and interpret the 
dyadic K [see Eq. (I14II ] that appears in the Lippmann- 
Schwinger equation for the field F: 



G^(r,r',c.) 



A 



£(r)(w/c)2 



^e(r',r). (21) 



It consists of the generalized transverse Green function 
(|20|l and a term proportional to the transpose of the gen- 
eralized transverse delta function, C3- Both terms 
are medium-dependent. Note that K is generalized trans- 
verse in its variable r. If only because of this property, 
K is not equal to the total Green function H16|l . Nev- 
ertheless, the definition Ijl9b|l of the longitudinal Green 
function can be used to rewrite K as 



K(r,r',c^) = G(r,r',c^) 



1 



e{Y){u/c)' 



r<5(r-r')l. (22) 



According to this identity, the dyadic K differs from the 
full Green function of the medium only when its two posi- 
tion arguments r and r' coincide. Although different from 
G, the quantity K will also be called a Green function. 
The occurrence of K rather than G in the Lippmann- 
Schwinger equation will be discussed further in the Ap- 
pendix, where the volume-integrated electric field around 
an atom is calculated. 



IV. SINGLE-ATOM PROPERTIES ALTERED 
BY THE MEDIUM 



A. Solution of the LS equation 

An atom in a group of atoms in an inhomogeneous 
dielectric will have different properties as compared to 
free space, because of the dielectric and because of the 
other atoms. In this section the effect of the medium on 
the individual atoms will be considered. The next and 
major step, in sectional will be to study some effects 
that the medium-modified atoms can have on each other. 

Assume that in the dielectric there is only one guest 
atom present with dipole moment /x and transition fre- 
quency O. The effect of the medium on the scattering and 
emission properties of the atom can be found by solving 



Eq. Hll() exactly by successive iterations 



e(") + k • s 



K V 



e(") + k • s 



K V K V E(°) + K S 



(23) 



In this equation, K and V are classical quantities, whereas 
F, E^*'', and S are quantum mechanical operators. The 
infinite series of multiple-scattering terms can be summed 
to give 

F(r,tj) =E(")(r,cj)+F,cat(r,cj)-KF,ourcc(r,tj), (24) 

where, as before, E(°^ (r, uj) is the electric-field operator 
of the inhomogeneous medium in the absence of the guest 
atoms. 

The operator Fscat(r, oj) in Eq. 124|) describes light that 
is scattered by the guest atom and it has the form 

F,cat(r, u) = K(r, R, u) ■ T(c^) • E^") (R, u), (25) 

with the single-atom T-matrix defined by 



T{lo) = jlT{uo)fi — jj, 



I- fi- K(R,R,tj) • fiV{uj) 



. (2^) 

The T-matrix is sometimes written as — (w/c)^ times a 
dynamical polarizability aioj) [compare with Eq. I|13|l ] 
and both depend on the atomic position inside the inho- 
mogeneous dielectric. The expectation value of the scat- 
tered field (|25l) only depends on the initial quantum state 
of the light (through the term E'^*'^). Unlike for a two- 
level atom, the light-scattering properties of a harmonic- 
oscillator atom do not depend on the atomic excitation. 
The scattering process can be read from right to left in 
the right-hand side of Eq. (^3)) : light E^^^ that has not yet 
seen the atom scatters off the atom (as described by T), 
and the scattered part of the light propagates through 
the dielectric as described by K. 

Finally, there is in Eq. H24|l the source field 

Fsourco(r,c^) = K(r,R,cj)-S(cj) (27) 
+ K(r, R, w) • J{uj) ■ K(R, R, uj) ■ S(w). 

Expectation values of the source field Fgourco only de- 
pend on the initial atomic state. Notice that the same 
T-matrix that shows up in the scattered field (|25|l also 
appears in the source field (|27|) . Light emitted by an 
atomic point source will be studied further in Sec. IIVBI 
and scattered hght in Sec. IIV Gl 



B. Light emitted by a point source 

The source field (|27|l can be rewritten as 
K(r,R,cj) • S(cj) 



3(r,a;) 



l-/i-K(R,R,w)-/iT/(w)' 



(28) 
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with S(w) as defined in Eq. H12() . The time dependence of 
the source field at the position r due to the presence of the 
source at R follows from the inverse Laplace transform 
of Eq. (EHJl, 



3(r,t) 



1 

2^ 



,(r,c^). (29) 



This integral can not be evaluated further without the 
explicit knowledge of the Green function K. The source 
field decays in time due to spontaneous emission by the 
atom. The decay rate can be found by multiplying nu- 
merator and denominator in Eq. (|28|l by {uj'^ — Vi^) and by 
realizing that the zeroes lo ~ of cj^ — — 2VlX{u}) 

with 



X{uj) = p, ■ K(R, R,Lu)-fi [^2wV(feoc^)] 



(30) 



are the frequency poles of the source field. Until now, our 
solution is exact. At this point we make a pole approx- 
imation, which is only valid if the atom-field coupling is 
weak. The pole approximation gives Qi = Q+X{fl), with 
X{^1) the difference between the dressed resonance fre- 
quency r2i(f2) and the bare atomic resonance frequency 
51. The exponential (amplitude) spontaneous-decay rate 
is 



r/2= -luiX{n). 



(31) 



The decay rate of the intensity of the field is F. It is non- 
negative by definition of K in Eq. (|14() . The delta func- 
tion term in K (122(1 and the longitudinal Green function 
G (R, R, n) in Eq. H19a|l are real quantities, so that F is 
proportional to the imaginary part of only the general- 
ized transverse Green function G""". The property that G'^ 
does not contribute to the spontaneous-emission rate is a 
generalization of the well-known result for homogeneous 
dielectrics ji^ and it only holds for non- absorbing di- 
electrics. Using the mode composition Eq. (|2UI) of G , we 
find T ~ n/ (feo) J2\ lA'-fA(r)|^WA'5(r2 — wa), the same ex- 
pression that one also finds from Fermi's golden rule 0. 
The decay rate depends both on the atom's position and 
on its orientation inside the inhomogeneous dielectric. 
For free space, the imaginary part of Gq (R, R, fl) is equal 
to — n/(67rc)l [see Eq. (|60af) ]. This gives the familiar free- 
space spontaneous-decay rate Fq — fj.'^V,^ / {3nheoc^). 

The dressed resonance frequency fii can be written as 
il + A'(ri) — ir{n)/2. Apart from a decay rate there is 
a frequency shift A' that is equal to Re X{n). For two 
reasons. A' is infinitely large even for free space. Firstly, 
the delta function term (5(r-R)l/[e(R)(w/c)2] inEq. ^ 
diverges when r and R are equal. This self-interaction 
term is medium-dependent through the factor e(R), but 
here and in the following we assume that guest atoms are 
electronically well separated from the dielectric medium, 
so that an empty-cavity model applies where the relative 
dielectric function is equal to unity at the position of 
the guest atom . The second reason why A' diverges 
is well known for free space: Go(r, R, w) diverges when 



r approaches R [see Eq. (|60|) ]. By a procedure called 
mass renormalization, the combined radiative shift in free 
space becomes finite, see for example [^-^j . From now on 
we can assume that 51 is the observable atomic frequency 
in free space; inside a dielectric, the atomic frequency 
shifts by an amount A that is given by the real part of 
[X{llj) — A"o(fl)], or in terms of the Green functions 



A = fi-Re[ G(R, R, fi) - Go(R, R, 51) ] • A 



heoc^ 



(32) 

The shift depends on the atomic position and dipole ori- 
entation. Notice that the full Green function is needed 
to determine the line shift, whereas for the decay rate it 
sufficed to know G""". 

The position-dependent radiative shifts are a mecha- 
nism of inhomogeneous broadening of the detected light. 
Electronic shifts usually dominate inhomogeneous broad- 
ening. Experimentally it will be hard to single out radia- 
tive shifts (a photonic effect) from electronic line shifts 
(due to changes in the atomic wave functions inside the 
medium) . 



C. Light scattered by a point scatterer 

In the scattered field of Eq. (|25|l , the atom appears as 
a point scatterer with an internal resonance in the opti- 
cal potential V{uj) and a corresponding resonance in the 
T- matrix in Eq. H2()|l . The scattered field has frequency 
poles in the T-matrix (just like the source field), but it 
also has poles for every optical mode frequency ±uj\ (un- 
like the source field). The time-dependence of the scat- 
tered field can be understood by separating the frequency 
poles (straightforward, but not spelled out here), again 
followed by an inverse Laplace transformation. In the 
following, place the atom in the origin. For the part of 
Fscat featuring the annihilation operators, one finds 



Fl+l(r,t) = E 



M /^A (0), 



ar(0)fA(O).Ax 



/oc 
dww2e-»"*K(r,0,w) - A X 
-OO 





251 




1 




+ A)2 -cj2 +F2/4- 


ILOxT LO 






51/(^^-H A) 


1 




sn 


-A — uj\^iT/2 uj - 


-51- A-f 


+r/2 




51/(51 + A) 


1 






-A + tJA + ir/2 uj^ 


-51 + A + 


+r/2 



(33) 



The negative- frequency part F^^^ of the field equals 
[F^+^Jt. The three terms between curly brackets in H33|) 
correspond to different optical processes. The first term 
describes elastic light scattering by the guest atom in- 
side the inhomogeneous dielectric; the second term has 
an exponentially decaying time dependence and corre- 
sponds to resonance fluorescence; finally, the third term is 
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an exponentially decaying nonresonant term, correspond- 
ing to an utterly improbable process that one could call 
anti-resonance fluorescence. In a rotating- wave approx- 
imation this process would disappear. After neglecting 
this third term, all the cj-poles in the integral H33|l have 
positive real parts, so that (|33(l can be called the positive- 
frequency part of the field Fgcat ■ 

Now consider the second term in Eq. (|33|l in more de- 
tail. In the resonance fluorescence process the guest atom 
is excited by light of frequency Wp, after which the atomic 
source decays exponentially due to light emission at fre- 
quency r^s. (In contrast, elastically scattered light oscil- 
lates with the pump frequency cUp.) The fluorescent light 
has the same position-dependent emission rates r(R, il) 
(|31|l and line shifts A(R, fl) (|32|l as found for sponta- 
neous emission before. A difference between the source- 
field of Eq. (jSHl and the fluorescent light in Eq. ((23) is 
that the latter also contains the information how well 
the pump light that comes in via mode X ~ p can excite 
the atom, in the factor fp(0) • /i. This difference is es- 
pecially important for inhomogeneous dielectrics, where 
atoms will be excited easier here than there. And indeed, 
it is through the process of resonance fluorescence that 
lifetimes of atoms in dielectric media are usually mea- 
sured. 

In a resonance fluorescence experiment, a light pulse 
or wave packet passes the atom during a time T. In ex- 
pression (|33|l . the intensity of resonantly emitted light 
depends on the expectation value with respect to the 
quantum state of light at time t — 0. This can only 
be a valid description of the process if T -C F^^, in other 
words, if the wave packet is so short that it "prepares per- 
cussionally the excited state" (^3> P- ^'^) of the atom at 
time t = 0. This is typically the case, even if the medium 
broadens the excitation pulse: excitation pulses last pi- 
coseconds and lifetimes lie in the nanosecond regime. 



V. SEVERAL ATOMS AS POINT SOURCES 
AND SCATTERERS 

A. Solution of the LS equation 

In section Irvl it was found how scattering by and emis- 
sion rates of single atoms are influenced by their dielectric 
surroundings. In the present section it is studied how the 
medium-modified atoms can influence each other. The 
atomic wave functions are assumed not to overlap each 
other and to be unaffected by the dielectric. The atomic 
positions can be arbitrary, so we can decide to choose 
the atoms on a line or on a lattice ^3 o^' at random 
positions. The general method to solve the Lippmann- 
Schwinger equation Ijlll) in this more complicated situa- 
tion is outlined here. In Sec. IVII the formalism will be 
used to study two-atom superradiancc inside an inhomo- 
geneous dielectric medium. 

For one atom, the Lippmann-Schwinger equation Hll|) 
was solved by summing a series to all orders of the atomic 



potential. In the present many-atom case all atomic tran- 
sition dipole moments, orientations, and frequencies are 
allowed to be different so that also all optical potentials 
Vj are different. The LS equation will now be solved by 
efficiently summing a somewhat more complicated series. 
Use the abbreviations F = r(r,a;), — F(Rm,Ci;), 
Km = K(r,Rm,w), Kmn = K(R,„, R„,aj), and intro- 
duce F^ = (r, a;) + ^„ K(r, R^, a;) • S^H. Also, 

F^^^ is shorthand for F(^^(R„). By iteration it follows 
that the field (F - F^^)) of Eq. ^ becomes 

N N 

EK V F^^) + K V K V F^^) 

■^n ' / ^ "^m *m "^mn »'n ^ ri 

n— 1 m,n— 1 

N 

-f K,„-V,„-K„p-Vp-Kp„-V„-F«-H(34a) 

m,p,n— 1 

This can conveniently be rewritten in terms of the single- 
atom T-matrices of Eq. (I26|l as 

N N 

Vk T f(i) + V K T K' T F^i) 

n—1 m.n—1 
N 

+ J2 K,„.T,„.K:„p.Tp.K;„.T„.F«+(.34b) 

m,p,n— 1 

Here the tensor is defined as (1 — Smn)K.mn, which 
by virtue of Eq. is equal to G'„„ = (1 — Smn)Gmn- 
A single-atom T-matrix already sums up all multiple 
potential-scattering off a single atom, which explains that 
neighboring T-matrices in terms of this series belong 
to different atoms. The equivalence of Eqs. (|34a|l and 
(j34bp can be seen by expanding single-atom T-matrices in 
terms of single-atom potentials. Now every higher-order 
term in Eq. Ij34b|l can be constructed from the previous- 
order term by inserting into the latter the N x N matrix 
with (z, j)-elements /i^ • G^^ - Tj ■ fij. By summing the geo- 
metric series of matrices and dropping the abbreviations, 
it follows that 

7{r,uj) = F(i)(r,tj) (35) 

N 

+ K(r,R„,c.).TW(a;).FW(R„,^), 

m,n— 1 

with the A^-atom T-matrix 

(36) 

The N X N matrix M(cj) is defined as 

Mij{uj) = [Sij - (1 - %)/ii • G(Ri,Rj,uj) ■ fLjTj{uj)] . 

(37) 

Eqs. (EHl) and ((23 neatly sum up infinitely many scatter- 
ing events which are not described by G. Light propaga- 
tion in between the scattering off one atom and the next 
one is described by G and need not be rectilinear, since 
G is the Green function of the inhomogeneous medium. 
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As before, the total field F consists of the part E^^^ 
that has not seen the atoms, a scattered part and a 
source- field part. The field operator that describes the 
scattering of light by the iV-atom system has the form 

N 

F,eat(r,c^)= K(r,R™,c^).T(^jH-E(°)(R„,c^). 

m.n—1 

(38) 

This is a generalization of the single-atom result of 
Eq. H25|l . It describes elastic scattering as well as res- 
onance fluorescence off N atoms in an inhomogeneous 
medium. The expectation value of Fscat depends on the 
initial quantum state of light only. Similarly, for the N- 
atom source-field that only depends on the initial atomic 
state, we find 

N 

Fsourcc(r,c^) = K''^nr,Rm,c^) • S™(^), (39) 



which generalizes Eq. H27|l . Here, K^^' is a Green func- 
tion of the inhomogeneous dielectric including the N 
atoms: 



KW(r,r',o.) = K(r,r', 



(40) 



N 



+ Y K(r,R„,c.)-T^j(c.)-K(R„,r',c.). 

m,n— 1 

For r and r' different from one of the positions R^, 
r',uj) is equal to what one would call the to- 
tal Green function G*-^' of both the dielectric and its 
guests. The solution (|39|) shows that the source-field 
K'-"'^-' (r, Rm, cj) • Sm(a;) that emanates from atom m is 
influenced by the positions, orientations, dipole moments 
and resonance frequencies of the (N-1) other atoms. No- 
tice that the same A^-atom T-matrix describes the N- 
atom source fields and scattered fields. The two-atom 
source field will be studied in Sec. IVII 



B. Interatomic interactions 

In the results of Sec. IV Al interatomic interactions can 
be identified. Before doing that, we briefly mention pos- 
sible interatomic interactions that we already neglected 
or that simply do not occur in our theory. In a minimal- 
coupling formalism there would be a direct atom-atom 
interaction in the Hamiltonian. In a multipole formal- 
ism, the only direct interaction between neutral atoms is 
an interatomic polarization energy |4^ . Classically, this 
interaction is zero unless the smallest spheres containing 
the atomic charges have nonzero overlap _42J. Quantum 
mechanically, this 'contact energy' is negligible unless the 
interatomic distance is of the order of the size of the 
atoms such that wave functions overlap. We assumed 
that the atoms were further apart. Together with the 
fact that atoms are much smaller than the wavelength of 



light, this allows us to make the dipole approximation in 
which atoms are considered as point dipoles. Therefore, 
direct interactions between the atoms are absent in the 
dipole Hamiltonian Eq. Our approximations make 
that the only interatomic interactions that we can find 
are retarded dipole-dipole interactions, mediated by the 
electromagnetic field. 

Indeed, in Sec. IV Al interatomic interactions showed 
up in the N-atom T-matrix as terms proportional to 
the causal Green tensor of the medium. For two non- 
coinciding positions Ri and R2, the interaction (with 
dimension: [frequency]) has the form 



J12 — «/(R-i,R2,^) 



^/xi • G(Ri,R2,w) • fi2 



heoc^ 



hso 



^/ii-fA(Ri)fi:(R2)-/i2^ 



(41) 



For the latter identity, Eqs. (|14|l and (|22|) were used. 
Only after making a pole approximation in Sec. lVII will it 
become fully clear why we identify precisely this expres- 
sion as the dipole-dipole interaction. Modes with eigen- 
frequencies lux = were absent in the dipole interaction 
(|lcf) and consequently are absent in the dipole-dipole in- 
teraction 1)41(1 . 

The Green function G can be written as the sum of 
the generalized transverse Green function G""" and a lon- 
gitudinal Green function G^ [recall Eqs. (HHJ-lEll)]- The 
dipole-dipole interaction can be split into two analogous 
parts. The generalized transverse part is 

7 fX} Ti \ MiM2t^^ Al • fA(R-i)fA(R-2) • /i2 

Jgtrans — ^— 



feo 



(42a) 

It is also called the 'resonant dipole-dipole interaction' 
(or RDDI) . The strongest contribution to this interaction 
comes from the modes A with eigenfrequencies ujx near 
uj, which explains the adjective 'resonant'. Notice that 
"/gtrans IS zcro whcu Lu is zero. The other part is the 
longitudinal dipole-dipole interaction Jiong that has the 
mode expansion 

•/long(Rl, R.2) = . f^(Rj)f*(R2) • p,2. 

net) ^ 

(42b) 

Notice that Jiong is independent of the frequency lo. It is 
the generalization of the static dipole-dipole interaction 
that is well known for free space. Both the generalized 
transverse and the longitudinal dipole-dipole interactions 
are given here in terms of generalized transverse modes. 
Both Jgtrans and Jiong are influenced by the medium. 

Both G"^ and G^ have nonretarded dipole terms, so 
that a change in a source term changes instantaneously 
the longitudinal and generalized transverse fields else- 
where. It is only their sum that is fully retarded. This is 
well known for free-space Green tensors (see ^^) and it 
holds likewise for the Green functions Gj^^j^ and Qr^a-m 
homogeneous dielectrics as given in the Appendix. 
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It might seem strange that the longitudinal interaction 
Eq. Ij42bp is not given in terms of longitudinal modes. 
The physical reason is that longitudinal modes do not 
couple to the atoms in our formalism (see Sec. Ill A|l . 
Still, apart from the generalized transverse solutions 
(with L0\ ^ 0) of the wave equation |5J), there are also 
longitudinal solutions q,y (with = 0). There is a 
mathematical identity that allows one to rewrite the lon- 
gitudinal interaction Eq. Ij42b|) in terms of longitudinal 
modes. The identity originates from the fact that the 
modes {f\, q^^} together span the entire space of functions 
h with / dre(r)|h(r)|^ < oo. This space consists of a 
subspace of generalized transverse functions and a longi- 
tudinal subspace. The completeness relation in the entire 

space reads Ea fA(i")fA(i"')£(i"') + E,. Qi'Wq^ (r')e(r') = 
(5(r — r')l, with I the unit tensor. It follows that for 
r ^ r', one can replace iii Eq. (|42b|l by 

minus the sum '^^(\^{y)c^^{t'). Incidentally, the longi- 
tudinal modes q,y of the medium are different from the 
free-space longitudinal modes, because of their different 
orthogonality relations / dre(r)qi,(r) • q*/(r) = 5^yi. 



VI. TWO-ATOM SUPERRADIANCE IN 
INHOMOGENEOUS MEDIUM 

The general results of section will now be applied 
to two identical atoms positioned in an inhomogeneous 
dielectric. Assume that the two atoms have identical 
electronic transition frequencies 17 and dipole moments 
/i = their dipole orientations (ii and (i2 need not be 
identical. The source field of this two-atom system is [see 
Eq. ^] 



Fsourcc(r,c^) = K^"^(r,Ri,c^).Si(^)+K(2)(r,R2,c^)-S2(c^). 

(43) 

The goal is now to calculate the Green function K*^^^ of 
the dielectric including the guest atoms, in terms of the 



properties of the medium and of the individual atoms. 

According to Eq. (|40|l . the Green function K*-^'' is 
known once the T matrix T''^^ (|36|I is determined; T^^^ 
can be found by inverting the 2x2 matrix M H37(l . in 
which the single-atom T-matrices occur that are given in 
Eq. (|26|1 and the Green function K of the dielectric in 
Eq. H22() . It follows that the two-atom T-matrix is 

-r(2) ^ ( /ii/ii/3Ti /tiAtzTiJiars 

1 - TiJf^T2/P^ V M2Mir2 JiaTi ^2/12/3^2 

(44) 

with the dipole-dipole interaction J12 defined in Eq. H41|) 
and (3 as ij.'^uj'^/ (feoc^). Each of the four matrix elements 
of T*-^-* is a dyadic of the same type as the single-atom 
T-matrix (|26|l . Now abbreviate K(r, Ri,a;) as K(rl) and 
similarly for other terms. The Green function K''^-'(rl) 
can be written with Eq. (|40|l as 



I 



-(2) 
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K(ll) 



r(2) 
12 



• K(21) 



K(^)(rl) = K(rl) 

+ K(r2) • [ Tg> ■ K(ll) + Tg^ ■ K(21) J . (45) 

Use Eq. H44() to rewrite the T-matrix elements of T*-^-* in 
terms of the single-atom T-matrices. The first one of the 
two parts of the source field (|43|l is associated with light 
initially residing in atom 1. This part can be written in 
terms of single- atom properties as 



K(^)(rl) ■ Si 



l + TiXi/l3 
l-T,J^,T2/f3\ 
K(rl) • Ai + K(r2) • fi2T2J2i/P ] m) 



with Xi — Xi{uj) as defined in Eq. (|3()|) . The source field 
has now been expressed in terms of the T-matrices of the 
individual atoms, but it is rewarding to break up the T- 
matrices in parts that depend on the medium alone and 
parts that depend on the atoms: 



K(2)(rl) - Si = 



(^2 _ n^)S-^ [ K(rl) • fii (^2 - f72 _ 2nX2) + K(r2) • fi2 2^.h2 ] 
(w2 _ f]2)2 _ 2n [Xi + X2) {uj^ - ^2) + 4172 1^X1X2 - JI2) 

I 



(47) 



The denominator carries the important information 
about the resonance frequencies Q.± (11) of the two-atom 
system. There are two resonance frequencies near Q, and 
two near — fi. If these resonance frequencies change lit- 
tle due to the electromagnetic coupling with the dielec- 
tric, we can replace the frequency dependent functions 
A'i,2(w) and Ji2(w) in the above expression by their val- 
ues in cj = ±17. This is the pole approximation that we 
also made for the single atom. We find the two resonance 



frequencies 

^^±(f^)=^^ + ^^±^(^^)'+^. (48) 

The other two resonance frequencies occur at — f7j_(f7), 
so that all four have negative imaginary parts. 

When the atoms are far apart, then J12 tends to zero 
and the two resonance frequencies 17± (|48|l are simply 
the two single-atom frequencies 17 1 and 172 with their 
medium-dependent radiative shifts Ai_2 and decay rates 
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ri_2- In the other extreme situation, for atoms with par- 
allel dipoles atoms (almost) on top of each other, — Imfi^ 
approaches twice the single-atom amplitude decay rate, 
whereas — ImJ7_ has the limiting value zero. Analogous 
to free space, $7+ corresponds to the superradiant state 
of the two-atom system in the medium, whereas f2_ is 
the frequency belonging to the subradiant state. 

Now rewrite Eq. (|47() as a sum over individual first- 
order frequency poles. With Eq. (|12|l one has 



K(2)(rl).Si ^ Y: 



[ K(rl) ■ fjLi[l± sin a) ± K(r2) • (12 cos a ] 



{uj + 9.)hi{Q) + {oj - VL)b\{Q) 
1 1 



(49) 



A (complex) angle a — Q;(f2) has been introduced which 
measures the inhomogeneity of the medium as felt by 
the two-atom system, in comparison with the atom-atom 
interaction: 



sma 



(r!i -f}2)/A, 



2J12/A, (50) 



with A equal to \J — 1^2)^ + 4J^2- When the angle a 
is zero (such as in free space), the atoms are said to be 
placed at equivalent positions in the medium. 

In the expression H48|l for the resonance frequencies and 
in the angle a (I5()|l . a driving term and a detuning can be 
discerned. The driving term is the dipole-dipole interac- 
tion J12 and it signifies how important the one atom is 
as a light source for the other. The term (ili — r22)/2 is a 
detuning: larger medium-induced local differences felt by 
the identical atoms make the resonant transfer of a pho- 
ton between them less probable. The driving term and 
the detuning have the same physical origin and can not 
be changed independently. By bringing the atoms much 
closer in each others near field, they will be tuned better 
and interact stronger at the same time. The outcome 
of the competition between medium-induced driving and 
detuning will be studied in an example in Sec. IVIII 

The time dependence of the source field can now be 
calculated with an inverse Laplace transformation. No- 
tice that the positive-frequency poles in Eq. 1491) have 
negligible residues in terms proportional to b\; similarly, 
negative-frequency poles hardly contribute to terms in- 
volving the annihilation operator 61 (0) and can be ne- 
glected as well. The total source field Fsourcc(r, i) is the 
field H49|) that originates from the initial excitation of the 
atom labelled 1, accompanied by the source field that 
originally came from the second atom: 

source 

(r, t) = Li(r, i)6i(0) + L2(r, i)&2(0) + H.c. (51) 
The vector Li can be written as the sum of Li_|_ and Li_ , 



with 

Li±(r,i) = 



dw 



X (52) 



[ K(r, Ri, w) • /ii(l ± sin a) ± K(r, R2, t^) • A2 cos a ] . 

This is the central result of this section. Vectors L2± can 
be found by interchanging the indices 1 and 2 in the right- 
hand side of Eq. H52() , which also causes a sign change in 
sina (f^ . 

Eq. (|52|l describes the full time dependence of the 
source that is excited at time i = 0. Initially light has 
been emitted but has not arrived at the detector yet. 
This initial phase lasts a certain delay time td, depending 
on the optical path length between source and detector. 
The initial phase is followed by a transient regime, in 
which light that has chosen the shortest path already ar- 
rives at the detector at r, while light that takes a longer 
path has not arrived yet. The transient regime can be 
neglected if it lasts much shorter than the typical atomic 
decay time, which is usually the case. Assuming the same 
delay time for both resonance frequencies and for both 
atoms, we find 

Lr±(r,t).::^^^Ml±M^(,_,,)e--.*x (53) 
4eoc^ 

[ K(r, Ri, r2) • /ii(l ± sina) ± K(r, R2, Vi) ■ (12 cosa ] . 

In this equation we see that the source amplitudes of 
the atoms are influenced by their environment. There 
are overall factors which are equal for both atoms. The 
atoms differ in that the source amplitude of the light that 
is finally emitted by the first atom has a factor (l±sina) 
while the corresponding factor for the second atom is 
cos a. The results for the source field can be inserted 
into the intensity operator 



/(r,i) = 2£ocE(-)(r,t) • E(+)(r,i), 



(54) 



to give the time-dependent intensity of the light emit- 
ted by the two atoms, at a detector position where e(r) 
equals unity. Suppose that the two atoms share a single 
excitation so that their initial state is the superposition 



*(t = 0)) = pb\{0) + e">'^\~p^hl{0) |0). (55) 



Then the expectation value of the intensity operator is 



(/(r,i)) = 2eoc^p'\U? + {I - v')\U\'' 

+ 2p^l - p2Re [ e'"*LJ • L2 ] |, (56) 

where variables (r, t) were dropped. 

From Eqs. (|53l) and H56I) it follows that both atoms act 
as sources and superradiance can take place, even if only 
the first atom is initially excited (so that p = 1). This 
point was stressed for two-atom emission in free space in 
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[49j. The time-dependent intensity that passes at r is a 
comphcated interference pattern of source fields emitted 
by four sources: a fast(er) and a (more) slowly decaying 
source at Ri, and also a fast and a slow source at the 
atomic position R2. The photon is shared by and ex- 
changed between the atoms until it is finally emitted, via 
either the fast superradiant or the slow subradiant decay 
process. Amplitudes of the fast and slow sources origi- 
nating from an initially unexcited atom depend both on 
the interaction between the atoms and on their medium- 
induced detuning. 

VII. APPLICATION: SUPERRADIANCE NEAR 
A PARTIALLY REFLECTING PLANE 

Our multiple-scattering formalism will now be applied 
to the situation of two identical atomic dipoles in the 
vicinity of an infinitely thin plane that partially reflects 
light. We are interested in medium-induced spontaneous- 
emission rates, Lamb shifts, interatomic interactions, and 
sub- and superradiant decay rates. These quantities of 
interest can be expressed in terms of the Green tensor of 
the medium. In a recent paper [H^], we already devel- 
oped a method to efficiently calculate the Green tensor 
of a medium consisting of one 'plane scatterer' and of 
several parallel planes, but we only used it to calculate 
spontaneous-emission rates. For a more detailed discus- 
sion of the plane-scatterer model and the method to cal- 
culate the Green tensor, we refer to Below, we give 
a brief outline of our calculations. We will then focus on 
those aspects of our results that we believe are generic 
for many more inhomogeneous dielectrics. 

The atomic dipoles are assumed identical and parallel 
to each other {fii = ^2), and parallel to the plane. More- 
over, assume the atomic positions — (xi,yi^Zi) to be 
confined to the line xi = iji = 0. The plane is assumed 
perpendicular to the z-axis. The interatomic interaction 
Eq. 141|) can then be written as 

J12 = J(Ri,R2, 1^) = To ^ dfc||A;|| (G^^ + G™) . 

(57) 

Here, G^'' stands for the component of the Green tensor 
G(fc|| , zi, Z2, ri) that describes propagation of s-polarized 
light, while G^" describes p-polarized light 5(J|. For a 
single plane, we have 

G^^(fc||,zi,z2,f^) =Gg^(fc||,zi,z2,f^) (58) 

-|-Gg^(fc|| , Zl, Zplane, 1^ )T^^ (fc|| , r!)Gg^ (fc|| , Zpianc , Z^,^). 

The free-space tensor component GQ*(fc|| , zi, Z2, f2) equals 
eyiY>{ikz\zi — Z2|)/(2iA:z), with the wave vector kz defined 

as ^(rj/c)2 - fc2. The T-matrix T^"(fc||,fi) of the plane 

for s-polarized light has the form — [(Deff (f^/c)^)~^ — 
i/{2kz)]~^ ■ The plane is fully characterized by the sin- 
gle parameter Deff: which we call its 'effective thickness'. 
We choose the value DcS = 0.23A. With this choice. 



32% of s-polarized light is transmitted through the plane 
when averaged over 47r incoming angles. Higher val- 
ues of -Deff give less transmission. For the Green ten- 
sor component G""" in Eq. (|57|) one can write an ex- 
pression analogous to Eq. (|58|1 : in the right-hand side 
of Eq. (|58|l . the components Gg" must be replaced by 
Gl'' = {kzc/nfG^Q", and T''^(fc||,rj) by r™(fc||,17) = 
-[(i^eff(f^/c)2)-i-zfc,cV(2r!2)]-i. 

If in the integral Eq. (|57|l the in-plane wave vector k\\ 
becomes larger than VL/c, then the wave vector kz be- 
comes purely imaginary and equal to in, with k equal to 

Y^fcy — (ll/c)2. The semi-infinite integration interval in 

Eq. (|57(l therefore falls apart into two parts: a radiative 
part with k\\ between and f2/c, and an evanescent part 
with fc|| from O/c onwards. The evanescent part of the 
integral is purely real, except that there is a purely imagi- 
nary contribution from a pole in T^* at k = I?eff (f^/c)^/2. 
This pole corresponds to an s-polarized guided mode. 
Near the pole, the real evanescent part of the integral 
over the s-wave integrand in Eq. (|57|l must be taken as a 
Cauchy principal-value integral. There is no correspond- 
ing pole in . The evanescent part of the integral for 
p-polarized light is purely real and can be evaluated nu- 
merically right away. All relative errors in our numerical 
results are smaller than 10~^. 

The single-atom spontaneous-decay rate Eq. H31|l can 
be found from the interaction Eq. H41(l through the 
relation r(Ri,fi) = -2ImJi2(Ri, Ri, f^), while the 
the Lamb shift Eq. follows from A(Ri,0) = 

Re[Ji2(Ri,Ri,l^) - j|2'(Ri,Ri,f^)]. Fig. HI shows how 
single-atom properties are modified by the presence of 
the plane. The figure shows a peak in the decay rate 
near the plane due to emission into the guided modes 
. Away from the plane, the decay rate shows damped 
oscillations towards the free-space decay rate Fq. There 
are two oscillations per wavelength A, a characteristic 
also well-known for spontaneous emission near a perfect 
mirror 44] . The Lamb shift shows similar damped oscil- 
lations around A = away from the plane. At distances 
less than A/10 the shift becomes strongly negative and it 
actually diverges to minus infinity. The atom is attracted 
to the plane |4J| , but here we assume atomic positions to 
be fixed. 

In Fig. 121 we present dipole-dipole interactions for two 
atoms near a plane. The first atom is kept fixed in the ori- 
gin, the distance of the plane to this first atom is chosen, 
and then the absolute value of the interatomic interac- 
tion Ji2 is plotted as a function of the position of the 
second atom, relative to the free-space value | \ . The 
interaction is the sum of radiative and evanescent inter- 
actions, of both s-polarized and p-polarized light. The 
figure shows that for Z2 approaching zi = 0, the relative 
difference between | J12I and the (divergent) free-space in- 
teraction strength | \ becomes negligible, irrespective 
of the position of the plane. This holds independently of 
the reflectivity of the plane (not shown in Fig.^l. Inter- 
estingly, the dipole-dipole interaction is also independent 
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FIG. 1; Spontaneous-emission rate F (solid line) and Lamb 
shift A (dashed line) as a function of the position of an atom 
near a partially reflecting plane. The effective thickness of 
the plane is Dctt — 0.23A. The atomic dipole moment points 
parallel to the plane. The plane is positioned at Zpianc = 0.4A, 
to make comparisons with later figures easier. Both F and A 
are given in units of Fo, and the positions z are scaled to the 
wavelength A = 2nc/fl of the emitted light. The period of the 
damped oscillations in both F and A is A/2. 



of the position of the plane (but not of its reflectivity) 
if the plane stands in between the two atoms. In other 
words, with the atomic positions fixed at either side of 
the plane, one can move the plane back and forth with- 
out changing the interatomic interaction. This fact can 
be read off from Fig. [21 for Z2/X > 0.4, where the three 
graphs (corresponding to three plane positions) overlap. 
It can also be understood from the form of the interaction 
in Eq. H57|l , because all terms in the interaction either de- 
pend on |z2 - zi| or on (|zi - Zpia„e| + \z2- Upland)- For 

Z2 < and \z2\ ^ A, the relative interaction |Ji2/'/i2''l 
approaches a constant value, which can be either larger 
or smaller than unity, depending on the distance of the 
plane to atom 1. As a check on our calculations (not 
shown), we found that interatomic interactions vanished 
(as expected) when an almost ideal mirror (a plane with 
Dcs = lOOA) is placed in between them. 

Fig.|3Ja) shows two-atom superradiant and subradiant 
decay rates, as modified by the presence of the plane. 
The plots are based on Eqs. (|^ and (|57l) . The com- 
plex square root in Eq. H48|l has solutions that differ by 
an overall minus sign. Care was taken to choose the 
solution from the same branch as we varied the posi- 
tion of the second atom. Without the plane, one would 
have had ro,± = Fo T 2Imj}2''. With the plane, the 
medium-induced detuning becomes negligible as Z2 ap- 
proaches zi. Then F_/Fo vanishes, as in free space. 
The corresponding small-distance limit of F_|_ /Fq is not 
equal to 2 as for free space, but rather twice the single- 
atom decay rate F = 1.14Fo in the presence of the plane. 



FIG. 2: Absolute values of interatomic interactions 
J12 (Ri , R2 , i^) near a partially reflecting plane, scaled to the 
free-space interaction strength | J^j' (Ri, R2, n)|. The plane 
is as in Fig.0 The first atom's position is fixed in the origin. 
The second atom travels along the line {x2 = y2 ~ 0,2:2). 
The atomic dipoles point in the same directions, parallel to 
the plane. The three graphs differ in the position of the plane 
with respect to the first atom. Solid line: Zpianc/A — 0.4; 
dashed line: Zpiane/A = 0.2; dotted line: Zpianc/A = 0.1. (All 
three planes are shown, but in each case considered only a 
single plane is present.) 



If the second atom moves towards the mirror, then the 
medium-induced detuning (see Fig. ^) grows fast while 
the dipole-dipole interaction (Fig. ^ decreases. With 
Eq. ||^ we then find that F„ ~ F2 and F4. ~ Fi. 
Indeed, for Z2 closer than A/ 10 to the plane, F_ fol- 
lows the single-plane emission rate of Fig. while F_|_ 
equals F(Ri,0) = 1.14Fo. With atom 2 so close to the 
plane, superradiance is completely absent, even though 
the identical atoms are less than half a wavelength apart. 
For Z2 > 0.5A or Z2 < 0.3A, detuning has become less im- 
portant and the decay rates follow (not quite sinusoidal) 
damped oscillations. Their period is A, as it is for super- 
radiance in free space. 

Fig. |3Ib) again shows super- and subradiant decay 
rates, now also for larger interatomic distances. For 
~4 < Z2/X < 3, the rates F± exhibit the same damped 
oscillations with period A that we also saw in fig. |3Ia) . 
However, for Z2/A smaller than -4 or larger than 3, F-t 
show two oscillations per wavelength, like we saw for the 
single-atom decay rate in Fig. ^ Hence we can identify a 
rather sharp cross-over regime at a few wavelenghts away 
from the plane between superradiance and single-atom 
emission. For larger distances, again medium- induced 
detuning dominates the dipole-dipole interaction. In- 
deed for large distances we see the same behavior as very 
close to the plane, namely that F_ approaches F2 (which 
at these positions almost equals Fq) while F+ has the 
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FIG. 3: Subradiant and superradiant decay rates for two 
atoms near a partially reflecting plane, as a function of the 
position of the second atom. The situation is as in Fig. |2| 
with the plane fixed at Zpianc/A = 0.4. Fig. (a) zooms in 
around the plane, showing a vanishing subradiant decay rate 
r_ (solid line) as the second atom approaches the first one in 
the origin. The superradiant decay rate r+ (dashed line) be- 
comes more than twice the single-atom decay rate Fq. Both 
F_ and F+ show perturbed oscillations on the scale of A. 
Fig. (b) zooms out to larger distances, showing a cross-over 
regime between damped oscillations with a period A and more 
distant oscillations with a period A/2. 



limiting value Fi — 1.14ro. In the cross-over regime, 
J12I has the same order of magnitude as the detuning 

lAi-A2-i(ri-r2)/2|. 

If one puts the plane closer to the first atom, then this 
atom becomes further detuned from its free-space prop- 
erties. The cross-over should then take place with the 
second atom at shorter distances where the interaction 
is still stronger. This we have verified (not shown). At a 
fixed frequency, the spatial intervals in which superradi- 
ance occurs therefore depend on three distances, namely 
the interatomic distance and the distances between each 
atom and the plane. 




FIG. 4: Absolute value of the source amplitude C2 = cos a 
of the second atom, divided by the source amplitude Ci+ = 
(1 -I- sina) of the first atom. The situation is as in Fig. |3 
The amplitudes are associated with superradiant emission, 
see Eq. (ICT . 



Not only the super- and the subradiant emission rates 
are influenced by the presence of the plane, but also the 
source amplitudes of the two atoms are modified, shown 
in Eq. (|53|1 : if initially only the first atom is excited, 
then the source amplitude of the second atom is modi- 
fied by a factor C\± = 1 ± sina, for superradiant (-I-) 
and subradiant (— ) decay, respectively. Atom 2 gets a 
factor C2 = cos a for both decay processes. Fig. ^ shows 
IC2/C1+I as Z2 is varied. The emission rates of Fig. Eland 
the relative amplitudes in Fig. 0] show the same cross- 
over regions between oscillations with periods A/2 and 
A. When the two atoms coincide, the detuning vanishes 
and IC2/C1+I equals unity. (For free space, IC2/C1+I 
equals unity everywhere, even if the atoms do not coin- 
cide.) The coinciding atoms are equivalent and superra- 
diance can occur. On the other hand, close to the plane 
at \z2/X — 0.4| < 0.1, atom 2 is strongly detuned and 
IC2/C1+I vanishes: the second atom emits none of the 
light initially residing in the first one and superradiance 
does not occur. 

Fig. 01 also shows that at larger distances (22 < —4 
or Z2 ^ 3), detuning is again strong enough to make 
emission by the second atom less probable than in free 
space. At these larger distances, the medium-induced 
detuning suppresses the net transfer of light from atom 
1 to atom 2 and superradiance does not occur. In the 
two intervals —4 < Z2A < 0.3 and 0.5 < Z2/\ < 3 
where superradiance does occur, we see that the peaks of 
the relative source amplitude IC2/C1+I are higher than 
unity. There the probability that light initially residing 
on the first atom is finally emitted by the second one is 
higher than in free space. The peaks of Fig. Q) corre- 
spond to positions of the second atom for which most 
light is finally emitted by the second atom, although ini- 
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tially only the first atom was excited. Interestingly, the 
peaks of IC2/C1+I become higher as the second atom 
moves away from the first. The highest peaks occur when 
the complex-valued dipole-dipole interaction (almost) ex- 
actly compensates the complex-valued detuning. (Such 
a resonant situation does not exist for identical atoms in 
free space.) For larger \z2\^ the dipole-dipole interaction 
becomes too weak to compensate for the detuning and 
Fig. 2] shows very abrupt transitions from superradiance 
to single-atom emission on both sides of the plane. 



VIII. CONCLUSIONS, DISCUSSION, AND 
OUTLOOK 

In this paper, a multiple-scattering theory was set up 
with at its heart the Lippmann-Schwinger equation Hll() 
that describes the electromagnetic field operators in an 
inhomogeneous dielectric with guest atoms present. We 
solved the LS equation exactly in terms of the properties 
of the atoms and the Green tensors of the medium, both 
when one and when several guest atoms are present. The 
solution for the electric field operator has three parts: 
a part that has not seen the guest atoms, a part that 
describes the scattering by the resonant atoms, and a 
part that describes the atoms as sources. 

Our formalism is a generalization of an already exist- 
ing point-scattering formalism for classical waves. The 
generalization is twofold: first, our formalism is valid not 
only for free space but for atoms in all dielectrics that 
can be described macroscopically in terms of a real rel- 
ative dielectric constant £(r). Second, it is a multiple- 
scattering theory in quantum optics rather than classical 
optics. In relation to this point we find the double na- 
ture of atoms both as scatterers and as sources of light. 
The formalism is quantum mechanical in the sense that 
it can describe the propagation and scattering of non- 
classical sources of light. These can be either external 
or atomic sources. In quantum optics, the medium must 
be described with more care, just like the quantum and 
classical descriptions of a beam splitter differ As 
for the beam splitter, classical light sources give classical 
measured fields in our formalism, since we described the 
guest atoms as harmonic oscillators. 

A nice feature of the LS equation 1)11(1 is that it fol- 
lows exactly from a dipole Hamiltonian that is the result 
of a canonical quantization theory. The Hamiltonian de- 
scribes guest atoms microscopically and treats the dielec- 
tric macroscopically. The atomic dipoles do not couple 
to the electric field operator E but rather to a field op- 
erator that we call F and that includes the atom's own 
polarization field. For free space this is a well known 
result. We find that the propagator for the field F in 
our LS equation is not the ordinary Green tensor G, but 
rather a Green tensor that we called K. There exists a 
simple relation (|22|l between G and K for an arbitrary 
dielectric. G can be split into the generalized transverse 
Green tensor G"^ that propagates the vector potential A, 



and the longitudinal Green tensor G . 

In the Appendix we showed that the volume-integrated 
electric field (|64|) produced in free space by an atomic 
dipole is equal to minus one third of its polarization field. 
This is an operator relation at finite frequency. A differ- 
ent (incorrect) relation would have resulted if the field F 
had been interpreted as the electric-field operator. We 
have not come across other work that addresses the re- 
lation between the dipole interaction, the occurrence of 
K rather than G in a multiple-scattering theory, and the 
volume-integrated electric field around a dipole. In this 
respect, our formalism also sheds new light on quantum 
optics in free space. 

The infinitely sharp single-atom resonance in the po- 
tential V obtains a radiative shift and a width in the T- 
matrix T. In our formalism, the position-dependent shift 
and decay rate are the summed effects of infinitely many 
light-scattering events off the atomic potential. The 
scattered-field operator for a single atom contains two 
parts: an elastic-scattering term and a term describing 
resonance fluorescence. Direct interatomic interactions 
are absent in the dipole Hamiltonian (fTc|l . Dipole-dipole 
interactions appear 'dynamically' in the solutions of the 
Lippmann-Schwinger equation for several atoms. An in- 
homogeneous medium modifies both the longitudinal and 
the generalized transverse dipole-dipole interactions, see 
Eq. gl- 

The multiple-scattering formalism has been used to 
study superradiance in an inhomogeneous medium. The 
often dominant electronic component to inhomogeneous 
broadening was neglected in order to focus on pho- 
tonic effects. As an application, we studied how dipole- 
dipole interactions and two-atom superradiance are influ- 
enced by a partially reflecting plane. We found position- 
dependent modifications of dipole-dipole interactions. 
For our choice of parameters, the plane suppresses su- 
perradiance if one of the atoms is very close or very far 
from the plane. Both atoms will then emit as if alone. 
For intermediate distances, two-atom sub- and superra- 
diance will occur. Due to the plane, emission rates are 
modified and so are the relative amplitudes of the atomic 
sources. Interestingly, we found that medium-induced 
complex detuning can lead to enhanced transfer of light 
from the one atom to the other, before superradiant emis- 
sion occurs. Also, we found sharp cross-overs between 
spatial intervals where superradiance occurs (with decay 
rates oscillating once per wavelength) and single-atom 
emission (two oscillations per wavelength). 

The length of the intervals in which superradiance oc- 
curs depends on the atomic positions with respect to each 
other and to the plane. This length could be called a 
"perpendicular coherence length" . This would comple- 
ment the concept of a transverse coherence length (or 
effective mode radius) [S^, |^. The latter concept is 
used in the analysis of cooperative emission in a planar 
microcavity when the atoms have the same z-coordinate, 
but have different coordinates in another direction. An 
important difference between the two lengths is that 
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only the perpendicular coherence length is influenced by 
medium-induced detuning. 

We believe that our results for cooperative emission 
near the plane are generic and that similar cross-over 
regions will occur in more complex dielectrics. Still, it 
would be interesting to study the influence of other di- 
electric structures on multi-atom processes, Bragg mir- 
rors for example, or 'optical corrals' |5^|53- Photonic 
crystals are also very interesting media, for which su- 
perradiance has only been studied in an isotropic model 
pH where all position dependence is neglected. 

Like near a plane, superradiance inside a real photonic 
crystal will be influenced by medium-induced detuning. 
As another application of our formalism, statistical dis- 
tributions of optical proximity resonances of many-atom 
systems can be studied, to find analogies and differences 
in inhomogeneous optical and electronic systems 54]. 

We made a pole approximation in a late stage of our 
formalism, after which we found exponential atomic de- 
cay. The pole approximation no longer holds when the 
atom-field interaction becomes strong [55| . The approxi- 
mation also breaks down if local densities of states jump 
steeply as a function of frequency near the atomic tran- 
sition frequency. There is a current debate whether pole 
approximations will break down at the band edges of 
realistic three-dimensional photonic crystals 56], like it 
is found for the isotropic model 15^^ . In principle, our 
formalism could also be used without making the pole 
approximation. 

Our theory is valid if frequency dispersion of the 
medium can be neglected. Now single- atom emission 
rates only depend on one frequency of the medium, so 
that dispersion is not important. On the other hand, ra- 
diative line shifts, interatomic interactions, and hence su- 
perradiant decay rates do depend on all frequencies of the 
medium. In our example of two atoms near a plane, the 
immediate vicinity of the atoms was free space. However, 
for atoms embedded in a e 7^ 1 part of a medium, line 
shifts would diverge unless frequency dispersion of the 
medium is taken into account |5a ]. This will also be the 
case for position-dependent radiative shifts in photonic 
crystals '5W\. It will be interesting to study the influence 
of frequency dispersion of the medium on cooperative 
atomic emission, for example based on Refs. [3lll60| . 
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APPENDIX: HOMOGENEOUS DIELECTRIC 

1. Delta and Green functions 

For a homogeneous dielectric with refractive index n, 
the property 'generalized transverse' reduces to trans- 
verse in the ordinary sense. The medium is transla- 
tional invariant, so that ^"^^'"(r, r') — ^'^''^(r — r'). The 
transverse and longitudinal delta functions appearing in 
Eqs. (|17|l and (jl9b|) now become 



i,,„,(r) = ^5(r)l--i^(l-3r®r) (59a) 
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'^homW = -5(r)l + ^(l-3r®r), (59b) 

where r is defined as r/]r], the unit vector in the direction 
of r. The sum of the transverse and the longitudinal delta 
function is simply S{r)\, since their 'dipole' parts cancel. 
Notice that n docs not enter these delta functions. The 
derivation follows the free-space treatment 48]. 

The dyadic Green function Ghom(r, r') = Ghom(r — r') 
for the homogeneous medium is the sum of a transverse 
and a longitudinal part. The transverse part is |4n| ] 



Ghom(r,t^) 



I - 3f (g) f 

47r(naj/c)^r'^ 

ainur / c 



(60a) 



[P{inujr /c)\ + Q{inujr/c)r f] 



with the function P{z) defined as (1 — z ^ + z ^) and 
Q{z) as (-1 + 3z"^ - 3z"2). With the use of the defi- 
nition H19a|l of the longitudinal Green function and the 
transverse delta function (|59all . the longitudinal Green 
function is found to be 



Ghom(l",w) 



Sir) 



I — 3f f 



Airinuj/cy^r^ 3{nuj/c)' 



-I. 



(60b) 



The delta-function term in GJ^^jj^ appears naturally and 
there was no need to add it 'by hand' as is done elsewhere 

MM. 



2. Volume-integrated dipole field 

The rigorous multiple-scattering formalism of Sec. ^ 
with the Green functions K will now be used to calcu- 
late the volume integral of the electric-field operator E in 
terms of the atomic polarization fields J2m -f*™ '^^ El' 
with the volume taken over a small sphere enclosing an 
atom. 

With the help of the Eqs. ((9H|) . and the defi- 

nitions of the source fields (|12|l and potentials (|13|) , the 
polarization field in frequency space can be related to 
other operators as 



Pm(w) 



) + V,„(tj) ■ F{Rm,Uj) ] 



(61) 
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There exists therefore a simple relationship between the 
field F and the polarization fields [use Eq. Hll|l ] 

2 ^ 

(62) 

This equation is still valid for all inhomogeneous di- 
electrics. Now assume that the sources are in free space. 
Consider the volume-integral of the field F over a small 
sphere (denoted by 0) containing only the source at R™, 
at its center. The integral is determined by the free-space 
Green function Ko(r, Rm,w) for positions r close to Rm 
(see Eqs. (|21f) and H22|) '). The transverse Green func- 
tion Gj(r — RmjO;) in Kq has a vanishing contribution 
to the integral, since its pole goes as |r — Rm| at short 
distances (see Ea. (|6Uall . The dipole part of the trans- 
verse delta function H59a() has a vanishing angle-integral 
over the sphere and does not contribute either. What 
remains is the delta-function part of the transverse delta 
function, which gives the radius-independent result 

/ drF(r,c^) = -^P™(c^). (63) 

Now the subtlety becomes important that the field F = 
— D(R)/[eo£(R-)] is equal to the electric field E every- 
where except at the positions of the guest atoms [see 
Eq. Q]. The expression in Eq. H63(l is therefore not 
equal to the volume-integrated electric field. With the 
definitions of the fields D and F given in Sec. Ill Al one 
obtains the relation 

/ drE(r,u;) = --^P„(a;). (64) 



The static and classical version of this 'sum rule' is pre- 
sented for example in [6ll |. There, and more recently in 
|4(tI | , a delta function is added by hand to the static dipole 
field or to the longitudinal Green function. In contrast, 
Eq. H64|l was found here as an operator relation without 
adding any terms by hand. 



The interpretation of the field to which a dipole cou- 
ples is not just a matter of choice in the present formal- 
ism. If one wrongly identifies F as the electric field E 
but correctly derives the relation 1211) or (|22|) between K 
and G, then the wrong volume- integrated electric field 
2Pm(a;)/(3eo) would have resulted. The delta function 
term that is the difference between G and K in H22|l and 
the difference between the field operators E and F have 
the same physical origin: the atomic polarization field. 



Still, there is nothing truly quantum mechanical about 
the sum rule H64|l . In a classical canonical theory, one 
would find the same dipole coupling —fi ■ F and Green 
function K. However, a canonical formalism is usually 
by-passed in classical optics. It is then assumed that a 
classical dipole couples to the classical electric field and 
furthermore that light propagates from a source accord- 
ing to the Green function G rather than K. By summing 
Eqs. H60a() and Ij60b|l for n — 1, one finds that Go na- 
turally has the correct delta-function term to produce 
Eq. 164|) . Therefore, although following a less rigorous 
procedure, one has the luck that there is no need to add 
terms by hand in order to derive Eq. 164|l classically. 
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